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Abstract 

The squashed Kaluza-Khen (KK) black holes differ from the Schwarzschild black holes with 
asymptotic flatness or the black strings even at energies for which the KK modes are not excited 
yet, so that squashed KK black holes open a window in higher dimensions. Another important 
feature is that the squashed KK black holes are apparently stable and, thereby, let us avoid 
the Gregory-Laflamme instability. In the present paper, the evolution of scalar and gravitational 
perturbations in time and frequency domains is considered for these squashed KK black holes. The 
scalar field perturbations are analyzed for general rotating squashed KK black holes. Gravitational 
perturbations for the so called zero mode are shown to be decayed for non-rotating black holes, 
in concordance with the stability of the squashed KK black holes. The correlation of quasinormal 
frequencies with the size of extra dimension is discussed. 



1 Introduction 



Recent years, higher dimensional black holes have become one of the key objects of the modern high 
energy physics. In brane world scenarios with large extra dimensions [I], it was suggested that mini 
black holes may be created in particle colliders [2j, which encouraged thereby detailed investigation of 
behavior of particles and fields around the black holes. 

Among important processes, the proper oscillations of higher dimensional black holes as a re- 
sponse to an external perturbations have been actively studied recently [3]-[l6]. These oscillations 
are governed by the so-called quasinormal modes, which have been in the focus of gravitational re- 
search during recent years because of a few reasons. First, quasinormal modes are expected to be 
an observed fingerprint of a black hole's gravitational waves with the help of a new generation of 
gravitational antennas. Second, for some asymptotically anti de Sitter space-times the quasinormal 
modes have interpretation in the dual conformal field theory (CFT) [17j, and help thereby to study 
some finite temperature field theory processes, such as hydrodynamics of quark-gluon plasma [18] 
in CFT. Finally, complete numerical investigation of quasinormal modes makes it possible to prove 
(in)stability for some black holes, when analytical proof is difficult [19\ [20l I21j . 

Usually, the asymptotic flatness of black hole spacetimes, which is an idealization of the boundary 
condition for isolated systems, is required in four dimensions. In higher dimensions, however, since the 
extra dimensions are expected to be compactified into a small scale, so-called Kaluza-Klein geometry, 
then the higher dimensional black holes have the asymptotic structure of the Kaluza-Klein type 0. It 
is interesting to study how the compactness of the extra dimension affects the quasinormal modes of 
the Kaluza-Klein black holes. 

The simplest example of five-dimensional black objects with the Kaluza-Klein geometry is the black 
string, the direct product of four-dimensional black hole and a circle. These objects look different 
from four-dimensional black hole only at sufficiently high energies, when Kaluza-Klein modes are 
excited |22^ [23| [2^ [25] . Therefore within these space-times we need high energy regime to see the 
extra dimension. 

If we allow the asymptotic structure of twisted S^ bundle over four-dimensional Minkowski space- 
time, there exist exact solutions of Kaluza-Klein black holes with squashed horizons for neutral |26ll27j . 
charged [28l [29], and extremely charged cases jSOj [31]. Such black holes look like five-dimensional 
squashed black holes near the event horizon, and like a Kaluza-Klein space-time at spatial infinity, i.e., 
locally a direct product of four-dimensional Minkowski space-time and a circle. Physical properties 
of the squashed Kaluza-Klein black holes are studied in [32\ [33l I34j . and generalizations of them 
appear in [35]. Owing to the non-trivial bundle structure, the size of the extra dimension might be 
observed even at low energies [33]. Then, it is interesting to study the the quasinormal modes of the 
squashed Kaluza-Klein black holes. This is even more motivated if we take into account that squashed 
Kaluza-Klein black holes are gravitationally stable [36j . 

In the present paper we study evolution of perturbations of squashed Kaluza-Klein black holes both 
in time and frequency domains and find the quasinormal spectrum for scalar field and gravitational 
perturbations of these black holes. 

The paper is organized as follows: in Sec II we study the quasinormal modes of scalar field 
perturbations around rotating Kaluza-Klein black holes with squashed horizon. Sec III is devoted to 
gravitational perturbations of non-rotating black holes. Finally in Sec. IV we discuss the obtained 
results. 

^ Asymptotically flatness would be a reasonable assumption for higher dimensional black holes if the size of horizon 
is sufficiently smaller than the size of extra dimensions. 
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2 Quasinormal modes of the scalar field for rotating squashed Kaluza- 
Klein black holes 



In this section, we investigate the quasinormal modes of the scalar field in the rotating uncharged 
squashed Kaluza-Klein black holes. 

The five-dimensional rotating squashed KK black hole with two equal angular momenta is de- 
scribed by[26l EH 
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r2 + a2 2 
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o"! = — sin tpdO + cos ip sin Odtf) , a2 = cos tpdO + sin ip sin 9d(/) , = dip + cos 
where < 6 < tt, < p < 2-it and < ip < in. The parameters are given by 
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Here n and a are parameters which correspond to mass and angular momenta, respectively, r = r+ 
and r = r_ are outer and inner horizons of the black hole and they relate to /x and o by = 
(r+r_)^,/i — 2a^ = + r^. The parameter Too corresponds to the spatial infinity. In the parameter 
space < r_ < r^- < Too , r is restricted within the range ^ r ^ ^oo* The shape of black hole horizon 
is deformed by k{r^). 

The wave equation for the massless scalar field ^{t,r,9, (pjip) in the background ([T]) obeys 



1 
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Taking the ansatz ^(t,r,6,(j),tp) 
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S{9), where S{9) is the so-called spheroidal 



harmonics, the radial and time variables can be decoupled from angular ones, so that the final wave- 
like equation reads 
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where / is the non-negative integer multipole number, \m\ < / and |2A| < 21 are integers, and 
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The radial coordinate p is given by 
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WKB (6th order) 


Frobenius 
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3.5149 


-0.15956i 


3.510564 - 


0.159674i 


10 


1/2 


3.5336 


- 0.15821f 


3.529225 - 


0.158317i 


10 


1 


3.5898 


- 0.15413f 


3.585417 - 


0.154231i 


10 


3/2 


3.6842 


- 0.14729i 


3.679778 - 


0.147368i 


10 


2 


3.8178 


- 0.13759i 


3.813421 - 


0.137647i 


10 


5/2 


3.9924 


-0.12493i 


3.988012 - 


0.124944i 


10 


3 


4.2103 


-0.10912i 


4.205922 - 


0.109076i 


10 


7/2 


4.4747 


- 0.08988i 


4.470523 - 


0.089745i 



Table 1: Scalar field perturbations. Comparison of 6th-order WKB and Frobenius approaches for the 
non-rotating uncharged squashed Kaluza-Klein black holes {po/p+ = 3, roo/r+ = 2). Quasinormal 
frequencies are measured in units of /9+. 



with 



2 ^0 ^ 2 , 2\ 



k{r = 0) = (10) 



»,4 
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Note that the three parameters po and p± = por'^/ (r^ — r^.) can define the metric (P) if Too < oo. In 
some papers they are used to parameterize the black hole instead of the parameters Too, r±. 
In frequency domain we used the following Frobenius expansion 



9 o s —IK oo / 9 9 \ 



~ ^ n=0 

where /t, v and $7 are chosen in order to eliminate the singularities at r = r+ and r = Too- The sign 
of K and is chosen in order to remain them in the same complex quadrant as io. We substitute 
(|lip into ([5]) and obtain the recurrence relation for the coefficients a„. After the recurrence relation 
is known we can find the equation with the infinite continued fraction with respect to which can 
be solved numerically |37j . 

The fundamental quasinormal modes will be of our primary interest hereafter, when dealing with 
frequency domain, because fundamental modes dominate in a late time oscillations. Contribution of 
all overtones can be seen in time domain considered in the next section. The fundamental modes for 
different values of A are shown in Fig(T]and Figl2l and tables 1 and 2. There one can see that the real 
oscillation frequency exerts some irregular growth (with local minimums) when Too is increasing until 
some moderately large values of r^. At larger r^o the growth of Re(w) changes into monotonic decay. 
The imaginary part of uj that determines the damping rate also has some initial irregular growth 
when Too increases, but at larger Too the two scenarios are possible: either monotonic decay (for large 
values of ratio r_/r+) or monotonic growth (for small and moderate r_/r+) (see Fig. 2). Thus, for 
a given mass and angular momentum of the black hole, one can learn the size of extra dimension 
from values of quasinormal modes of the emitted radiation. 

The calculations performed with the help of the Frobenius method are accurate and WKB values 
given here are rather an additional check. From tables 1 and 2 we can see that for small values of 
ratio Po / P+ the WKB method extended until the 6th WKB order beyond the eikonal approximation 
|381 [39] is in a very good agreement with the Frobenius method, but the larger po/ p+ the worse 
convergence of the WKB series. 
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Figure 1: Real parts of quasinormal frequencies as a function of r_ and r^o (/ = 0). 




Figure 2: Imaginary parts of quasinormal frequencies as a function of r_ and r^o {1 = 0). 



5 



3 Gravitational quasinormal modes for non-rotating squashed Kaluza- 
Klein black holes 



In the previous section, we have considered the scalar field. However, the tensor perturbations are 
more interesting from the point of view of the stability and the gravitational waves from black holes. 

The metric of the uncharged non-rotating squashed Kaluza-Klein black hole is a particular case 
of the previous metric 

ds' = -F{p)dr' + ^dp' + + (12) 

where we have difined r = 2pQt/rao and 

F{p) = l-^, G{pf = l + ^, rl = 4po{p+ + po). (13) 



Here, we have used a basis 



a^ = ha'Tia'). (14) 



The perturbed metric is 

h = hABdx^dx^ + 2hAidx^a^ + hija'a^ (15) 

where indices run A = t,p, i,j = +,—,3. Since the spacetime has the symmetry SU{2) x U{1), 
the metric perturbations can be classified by eigenvalues J,M for SU{2) and K for U{1). Here we 
consider only zero modes J = M = 0. Even in this case, since carry eigenvalues K = ±1, each 
components could have different eigenvalue K. It is important to recognize that the components with 
different K are decoupled. That is why we have the master equation for each K. To obtain master 
equations, we choose the gauge condition as 

^3+ = hs- = = htt = ht3 = (16) 

As is shown in [36j, the perturbation equations for the |iir| = 2 mode can be reduced to the wave 
equation for with the effective potential in the form 

jp - P+) 

Wp^po{p+ + po){p + po)^ 

+8p^Poi2pl - 5p+po + 25pl) + ppli20pl - 9p+po + 35pl) + 9p+plip+ + po)] . (17) 



6V + 256/5^0 - 32p>o(/5+ - H/^o) 



Similary, the perturbation equations for the = 1 mode can be reduced to the wave equation for 
/ip+ with the effective potential 

Vi = 3 , .3, ^I'^'w 2^, -,[i6p' + 128p'po-27pIpI{p++po) 

16p'^Po{p + Po) {p+ + po){p'' + 2ppo- p+po) L 

+32/po (P+ + 15po) + 2pV+/Oo (111P+ - 53/3+Po - 228p^) 

+8/po (-16p^ - 9p+po + IWpo) + PpIpI (-20/9^ + 151p+po + 187p§) 

Vpo (-827p^ - 959p+po + 764p2) + p^ p^ (-484p^ - 501p+/)o + 1103^^) 

+2p3p3 (36^^ _ miplpo - mp+pl + 126p^) 1 . (18) 



And, for the K = Q mode, we obtain the wave equation for /i33 with the effective potential 

{P - P+) 

16p3(p + po)3(4p + 3po)2 

+60p(13p+ + 12p)p[] + 9(9p+ + 35p)pg . (19) 



256/3+p^ + 64p3(l7p+ + 2p)pQ + 48p2(32p+ + llp)p^ 
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Po 


UJ 





UJ = 


0.226 - 


0.1711i 


0.1 


UJ = 


0.234 - 


0.2046i 


0.2 


UJ = 


0.243 - 


0.2003 


0.5 


UJ = 


0.263 - 


0.1848i 


1 


UJ = 


0.274 - 


0.1635i 


2 


UJ = 


0.271 - 


0.1344i 



Table 2: Perturbations of the K = type. Quasinomal frequencies for small p^. 



We shall use the time-domain analysis based on the standard integration scheme for the wave-like 
equations 

"^'^^""^ = -V^ip)Mt,p) , (i = 1,2,3) , (20) 

which is described for instance in [5l]. In detail, we applied a numerical characteristic integration 
scheme, that uses the light-cone variables u = t — and v = t + p^. In the characteristic initial value 
problem, initial data are specified on the two null surfaces u = uq and v = vq. The discretization 
scheme we used, is 

HN) = nW) + HE) - nS) - ^^nWMW) + ViEME) ^ ^^^^ 

8 

where we have used the following definitions for the points: N = (u + A,v + A), W = {u + A,ij), 
E = {u,v + A) and S = {u,v). We also used in this paper the WKB method developed in [35] and 
extended to the sixth order in [39]. As it is described in too many papers recently we refer a reader 
to the above mentioned papers for a detailed description. 

The quasinormal frequencies were obtained by WKB formula (Fig. 3) and by the above time 
domain integration (Table 3.). There are two distinctive features of gravitational perturbations, 
on the contrary to test scalar field perturbations. First, the K = 1 potential contains the small 
negative gap near the event horizon, which nevertheless does not induce any instability [36]. Second, 
K = 1,2 perturbations for some values of po have effective potentials without maximum but with 
a monotonic growth instead. For the latter case the WKB method cannot be applied but the time 
domain integration can be. 

From Fig. 4, 5 and 6, one can see that the quasinormal ringing for K = case and for K = 1,2 
cases are completely different. For K = perturbations one can see the usual damped oscillations 
until sufficiently large time, when they are dominated by the asymptotic late-time tails (see Fig. 4). 
For K = 1,2 cases we see that "tail behavior" occurs at much earlier time (see Fig. 5 and 6), as it 
takes place for perturbations of massive scalar field [42] . It is not surprise if one remembers that the 
effective potential for the massive scalar field with a large mass also has monotonic behavior as in our 
cases for = 1,2 gravitational perturbations. 

For larger values of po both real and imaginary parts of uj decrease (see Table 3). We can see also 
that for not large values of po the quasinormal modes of the Kaluza-Klein balck holes are longer lived 
(see Tables 3). 

4 Discussion 

In the present paper we considered the quasinormal spectrum of the scalar and gravitational per- 
turbations for squashed Kaluza-Klein black holes. The quasinormal frequencies of the scalar field 
are different from those of the Schwarzschild black holes and have smaller real oscillation frequencies 
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Figure 3: Real and imaginary parts of the gravitational fundamental quasnormal modes (n = 0) for 
K = (blue), K = 1 (red), K = 2 (yellow) modes. 



PO 


K = 


(t-d) 


K = (WKB) 


K = 


1 (t-d) 


K = 1 (WKB) 


K = 


2 (t-d) 


K = 2 (WKB) 


4 


0.2486 


- 0.1023i 


0.2467- 0.1019i 


0.332 


- 0.069i 


0.323 


- 0.062i 


0.506 


- 0.026i 


0.505- 


0.028i 


6 


0.2259 


- 0.0840i 


0.2255- 0.0842i 


0.282 


- 0.068i 


0.269 


- omoi 


0.380 


- 0.031i 


0.381- 


0.033i 


8 


0.2080 


- 0.0723i 


0.2083- 0.0728i 


0.249 


- 0.064i 


0.236 


- omii 


0.313 


- 0.034i 


0.315- 


0.035i 


10 


0.1943 


- 0.0644i 


0.1944- 0.0647i 


0.226 


- 0.058i 


0.215 


- omsi 


0.272 


- 0.035i 


0.273- 


0.035i 


12 


0.1828 


- 0.0585i 


0.1828- 0.0586i 


0.210 


- 0.051i 


0.198 


- omoi 


0.240 


- 0.034i 


0.244- 


0.034i 


15 


0.1686 


-0.0518i 


0.1688- 0.0519i 


0.189 


- 0.050i 


0.180 


- 0.046i 


0.212 


- 0.033i 


0.214- 


0.033i 


20 


0.1506 


- 0.0436i 


0.1511 - 0.0443i 


0.166 


- 0.045Z 


0.158 


- 0.040i 


0.181 


- 0.027i 


0.181- 


0.031i 


25 


0.1382 


- 0.0381i 


0.1380- 0.0391i 


0.150 


- 0.039i 


0.143 


- 0.036i 


0.160 


- 0.024i 


0.160- 


0.030i 


40 


0.1143 


- 0.0299i 


0.1127- 0.0302i 


0.120 


- 0.032i 


0.115 


- 0.029i 


0.122 


- 0.021i 


0.124- 


0.025i 


60 


0.0953 


-0.0241i 


0.0938 - 0.0242i 


0.100 


- 0.025i 


0.095 


- 0.023i 


0.098 


- O.OlSi 


0.100- 


0.021i 



Table 3: Quasnormal frequencies of K=l and K=2 gravitational perturbations of the non-rotating 
uncharged squased Kaluza-Klein black holes. Frequencies are measured in units of p+. 
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Figure 4: Time domain profile for gravitational perturbations = 0, /?o = p+ = 1- 
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Figure 5: Time domain profile for gravitational perturbations K = 1, pQ = = 1. 
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Figure 6: Time domain profile for gravitational perturbations K = 2, pQ = p-^- = 1. 
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and longer lived for not very larger pQ. The obtained evolution of gravitational perturbations shows 
damped oscillations in time domain, so that no instability is observed in concordance with analytically 
proved stability in [36]. Let us note that we analyzed here only "zero mode" perturbations, so that, 
strictly speaking, the opportunity for instability is remained in higher multipole perturbations. Yet 
we do not expect instability at higher multipoles because, unlike K = 1 mode, K = 2 mode does 
not have negative gap, so that higher K seem simply to increase the height of the potential barrier 
and stabilize the system. It is interesting that if we know the quasinormal frequency we can find the 
size of the extra dimension in the considered black hole model, so that quasinormal modes give a 
kind of opportunity to "look into" an extra dimension at low energies. In detail, when a dominant 
quasinormal mode is measured, then one can look into tables and plot numerical data in this paper 
and find out which is the value of po and the radius of the event horizon that corresponds to the 
observed quasinormal mode. In this way we can determine the parameters of the black hole and the 
size of the extra dimension, assuming that there exists no other Kaluza-Klein black holes with similar 
features. 

Our next step in the line of this research is to consider gravitational perturbations of the rotating 
squashed Kaluza-Klein black holes. Charged Kaluza-Klein black holes, especially maximally charged 
cases, are also interesting targets for investigating quasinormal modes, because the quasinormal fre- 
quencies of the Maxwell field and the gravitational field coincide in four-dimensional super-symmetric 
black holes [l3]. To check the stability of these black holes is an important problem. As the potential 
is frequency dependent and cumbersome, one definitely needs quasinormal modes analysis to test 
stability. 
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